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Abstract 

A "fat slit" is a compact domain in the upper half plane bounded by a curve 
with endpoints on the real axis and a segment of the real axis between them. 
We consider conformal maps of the upper half plane to the exterior of a fat slit 
parameterized by harmonic moments of the latter and show that they obey an 
infinite set of Lax equations for the dispersionless KP hierarchy. Deformation 
of a fat slit under changing a particular harmonic moment can be treated as a 
growth process similar to the Laplacian growth of domains in the whole plane. This 
construction extends the well known link between solutions to the dispersionless 
KP hierarchy and conformal maps of slit domains in the upper half plane and 
provides a new, large family of solutions. 
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1 Introduction 



Parametric families of conformal maps in 2D are known to be closely related to long wave 
limits of nonlinear integrable PDE's and their infinite hierarchies. This observation was 
first made in [1] for mappings of slit domains and then extended to mappings of domains 
bounded by Jordan curves in [2J. In both cases conformal maps from a standard reference 
domain (such as upper half plane, or unit disk) to a domain of a varying shape serve as 
Lax functions of an integrable hierarchy whose flows are identified with variations of the 
conformal maps described by an infinite set of Lax equations. The integrable structures 
arising in this way are dispersionless Kadomtsev-Petviashvili (dKP) and dispersionless 
2D Toda (dToda) hierarchies and, more generally, the universal Whitham hierarchy first 
introduced P H] in an absolutely different context with the aim to describe slow modu- 
lations of exact solutions to soliton equations. 

The most promising progress along these lines was achieved in geometric and physical 
interpretation of the dToda hierarchy. The key fact, established in |5] and further elabo- 
rated in [6l [7] , is that variations of domains under the Toda flows go exactly according to 
the Darcy law specific for growth processes of Laplacian type and viscous hydrodynamics 
in the Hele-Shaw cell with zero surface tension (see, e.g., [8l[9]). 

Although the dKP hierarchy is simpler than the dToda hierarchy, its role in the 
theory of conformal maps and Laplacian growth is not well understood. The connection 
with conformal maps observed in [T] gives a geometric interpretation to only rather 
special (degenerate) solutions of the dKP hierarchy which are reductions to systems of 
hydrodinamic type with a finite number of degrees of freedom. As is shown in [1] (see also 
plot [TT]). they are related to conformal maps of the upper half plane with slits emanating 
from the real axis. 

The aim of this paper is to demonstrate that the geometric interpretation of the dKP 
hierarchy is not limited by domains of such special kind. We show that the same dKP 
hierarchy is able to cover a much broader class of domains which can be obtained from 
the upper half plane by removing not just an infinitely thin slit but a whole compact 
piece (of a non-zero area and arbitrary shape) attached to the real axis, which we call 
a "fat slit" to stress the analogy (Fig. 1). There is an important difference, however. 
The dKP-evolution of usual slits is actually finite dimensional because the slits are to be 
regarded as arcs of fixed curves [H |TD|, so only their endpoints can move. In contrast, 
the dKP-evolution of fat slits (given by the Lax equations) takes place in an infinite 
dimensional variety corresponding to changing their shape in an arbitrary way. 

Let us recall the Lax formulation of the dKP hierarchy. Starting from a Laurent series 

oo 

zip)=p + Y.akp-' (1.1) 

k=l 

one introduces the dependence on an infinite number of "times" Ti,T2,T3, . . . via Lax 
equations 

dzip) . , . , . dBk dz dBkdz 

where the generators of the flows are polynomials in p of the form Bk{p) = {z^{p))^Q 
(polynomial parts of z'^(p)). The dKP hierarchy is an infinite system of nonlinear PDE's 
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Figure 1: a) A sUt; h) A "fat sUt" . 



for ttj's resulting from comparing coefficients in front of different powers of p in the Lax 
equations or in an equivalent system of equations of the Zakharov-Shabat type 



for all j,k > 1 (the Poisson bracket is defined in (11.21) ). The coefficients a, and the times 
Ti are assumed to be real numbers. 

Assuming that z{p) is a normalized conformal map from the upper half plane onto 
the exterior of a fat slit, we show that it obeys equations (11.21) . where Tk are properly 
defined harmonic moments of the fat slit (or rather of its exterior). This fact follows 
from the Hadamard formula for variations of the Green function with Dirichlet boundary 
conditions. In this sense the arguments are parallel to [T^l US] . Evolution in T = Ti with 
all other times fixed has an interpretation as a version of the Laplacian growth in the 
upper half plane with fixed real axis. 

Applying the approach developed in [121 1131 E] to the case of fat slits in the upper 
half plane, we construct the dispersionless "tau-function" (which is actually a limit of 
properly rescaled logarithm of a dispersionfull tau-function) of the dKP hierarchy as a 
functional on the space of fat slits explicitly given by 



It has a clear electrostatic interpretation as Coulomb energy of a fat slit filled with 
electric charge of a uniform density in the presence of an infinite grounded conductor 
placed along the real axis. This functional regarded as a function of harmonic moments 
obeys a dispersionless version of the Hirota relation which serves as a master equation 
generating the whole dKP hierarchy. 

2 Fat slit domains, their conformal maps and Green's 
functions 

Consider a compact simply connected domain B in the upper half-plane H bounded by a 
non-self- intersecting analytic curve 7 in H with endpoints x_, x+ on the real axis and a 
segment of the real axis between them. This segment will be called the base of B. Without 
loss of generality, one can assume that the origin belongs to the base. For brevity, and in 



+ {B,{p),Bk{p)} = 



(1.3) 




fat slit 



(1.4) 
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Figure 2: A fat slit B in the upper half plane and the corresponding symmetric domain 
D = BU B. 



order to emphasize an analogy with sht domains manifested in the common integrable 
structure of their conformal maps, we call such a domain a fat slit. Accordingly, the 
complement, EI \ B, will be referred to as a domain with a fat slit, or simply a fat slit 
domain (in our case, the fat slit half-plane). 

It is often convenient to treat fat slits as upper halves of domains symmetric with 
respect to the real axis. Namely, set D = B U B, where B is the domain in the lower half 
plane which is obtained from B by complex conjugation z ^ z (Fig. [2]). Obviously, the 
domain D is symmetric with respect to the complex conjugation. In what follows we call 
such domains simply symmetric. Vice versa, any compact simply connected symmetric 
domain D is a union of a fat slit and its complex conjugate. The boundary of D is assumed 
to be analytic everywhere except the two points on the real axis which are allowed to be 
corner points. 

2.1 Conformal maps 

Let p{z) be a conformal map from H \ B (in the z-plane) onto H (in the p-plane) shown 
schematically in Fig. [31 We normalize it by the condition that the expansion of p{z) in 
a Laurent series at infinity is of the form 



(a "hydrodynamic" normalization). Assuming this normalization, the map is unique. 
The upper part of the boundary, 7, is mapped to a segment of the real axis 
while the rays of the real axis outside B are mapped to the real rays [—00, p-] and 
00] (Fig. [3]). From this it follows that the coefficients Uk are all real numbers. The 
first coefficient, Ui := u, is called a capacity of B. It is known to be positive. We also 




fc>2 
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Figure 3: The conformal mapp{z). 



need the inverse map, z{p), which can be expanded into the inverse Laurent series 

u °° 

z{p) =p-- + J2 akP''' , bl ^ oo (2.2) 

with real coefficients connected with by polynomial relations. The series converges 
for large enough \p\. 

According to the Schwarz symmetry principle, the function z{p) admits an analytic 
continuation to the lower half plane. This analytically continued function performs a 
conformal map from the whole complex plane with a finite cut on the real axis from p_ 
to p+ onto the exterior of the symmetric domain D = B U B . 



2.2 Green's functions 

Let D = BuB be a symmetric domain and let G{z^ z') be the standard Green's function of 
the Dirichlet boundary problem in C\ D. The function is harmonic in C \ D with respect 
to both variables except aX z = z\ where it has a logarithmic singularity Gi^z^ z') = 
log \ z — z'\ -\- . . . and equals zero when either z or z' lies on the boundary of D. The 
Green's function solves the Dirichlet boundary value problem: the formula 

f{z) = -^i f{i)dnG{z,0\di\ (2.3) 

harmonically extends the function /(^) from the contour to its exterior. Here and 
below, dn is the normal derivative at the boundary, with the normal vector being directed 
to the exterior of D. We also note the Hadamard formula [15j for variation of the Green's 
function under variation of the domain: 

(5G(a, b) = — (f dnG{a, z)dnG{b, z)6n{z)\dz\ , (2.4) 
2n JdD 

where 6n{z) is the infinitesimal normal displacement of the contour. Some care is needed 
to define a deformation near the corner points. However, for our purposes it is enough 
to consider deformations with fixed corner points, then 6n{z) is well defined at any point 
of the boundary. Also, in this paper we consider only the case when both angles a+ and 
a_ are acute, < a± < 7r/2, then the normal derivative of the Green's function vanishes 
at the corners and the integral converges (see more details below). 

For symmetric domains the Green's function obeys the property G{z,z') = G{z,z'). 
Let us call a function f{z) even (respectively, odd) if f{z) = f{z) (respectively, f{z) = 
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—f{z)). It is natural to introduce even and odd Green's functions such that z') = 

G^{z,z') = G{z,z')±G{z,z'). 

Note that G~ {x, z) = for real x. The Poisson formula (12.31) for even and odd boundary 
functions can be written in the form 

fiz) = -^ I f{i)d^G^{z,i)m, (2.5) 

where the integration goes over the non-c/osec? contour 7 (Fig. Hj). The Hadamard formula 
for G^, 

5G^{a,h) = — I dnG^{a,z)dnG^{h,z)5n{z)\dz\ (2.6) 

directly follows from (12.41) . taking into account that deformations of symmetric domains 
obey the condition 6n{z) = 6n{z). 



2.3 The odd Green's function 



An important part in what follows is played by the odd Green's function G^ . It solves the 
following Dirichlet boundary value problem in H: To find a harmonic function in HI \ B 
bounded at infinity such that it is equal to a given function on 7 and on the rays of 
the real axis outside B. Similar to the Green's function G, G^ can be expressed through 
a conformal map to a fixed reference domain. The most natural reference domain in our 
case is the upper half plane H. It is easy to see that 



G {z, z!) = log 



p{z) -p{z') 



p{z) -p{z') 



(2.7) 



where p{z) is the conformal map (12. ip from H \ B onto H. We also need a useful formula 
for the kernel dnG~{a, z) in (12.51) through the conformal map, 



dnG [a, z) = —2Xmp{a) 



z E J 



\p{z) -p(a)|2 ' 

(which straightforwardly follows from (12. 7p ) and its limiting case as \a\ —>■ 00: 



dnlmp{z) = \p'{z) 



z e 'J . 



(2. 



(2.9) 



Let us present the expansion of the odd Green's function G (a, z) as \a\ 00: 

G-{a,z)=2j2llm{a-')Im{Bkip{z))). (2.10) 
k>i 

Here Bk{p) are Faber polynomials of p{z) defined by the expansion 



log 



p{z) -p 



E -1-Bk{p) 



2; — s> 00 



k>l 



and explicitly given by 



Bu{p) = {z\p))>o 



(2.11) 



(2.12) 
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where (. . .)>o means the polynomial part of the Laurent series. Indeed, fixing a point 
G EI \ B, we have 



>0 



log 1 - 



z{p) 



k>i k>i kz^{pi 

where pi = p{zi). To separate non-negative part, we write 

z{p) \ 



z{p 



iVJ>o 



log 1 



Z{j)i 



1 Pi - P , 1 z{pi) - z{p) 
log — — - + log ■ 



z{pi 



Pl-p 



and notice that the expansion of the first (second) term contains only non-negative (re- 
spectively, negative) powers of p. Therefore, 

^ ^ ^k{p) = - log(p(^i) -p)+ log^i, 

A:>1 

which coincides with (12. lip . In particular, Bi{p) = p. Clearly, 



BMz)) = z'' + 0{z-'), 
and the function Bk{p{z)) — z^ is analytic in C \ D. 



\z\ oo , 



(2.13) 



3 Local coordinates in the space of fat slits 

We are going to show that the harmonic moments Tk (defined as in (13. ip below) locally 
characterize a fat slit in the following sense. First, any small deformation of B that 
preserves the moments T^, is trivial, i.e., any non-trivial deformation changes at least one 
of them. This fact means local uniqueness of a fat slit having given moments. Second, 
the moments T^, under certain conditions discussed below, are independent quantities 
meaning that one can explicitly define infinitesimal deformations of B that change any 
one of them keeping all other fixed. In this weak sense they serve as local coordinates in 
the space of fat slits (c.f. the remark in Section 2.1 in [S]). 

3.1 Harmonic moments 

Given a fat slit B, let us introduce harmonic moments of the fat slit domain H \ B as 

2 r 

Tk = —Im / z-^d^z , k>2, 
nk JH\B 

(3.1) 

Ti = Im / z~^(fz . 

71 Jb 

Here we assume that the base of B is a segment containing zero. Since Im{z~^) < for 
z G H, Ti is always positive. Although the integrand in the formula for Ti is singular at 
the origin, the integral converges. The integral for T2 diverges at infinity, so one should 
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Figure 4: The integration contour 7. 



introduce a cut-off at some large radius and make the angular integration first; this 
prescription is equivalent to the contour integral representation given below. Note that 
this set of moments does not include the area of B. Note also that the standard harmonic 
moments dealt with in [13] are, for symmetric domains, real parts of the integrals in (13.11) 
rather than imaginary ones. 

Some other integral representations of the moments (13.11) are also useful. Imaginary 
parts of the integrals can be taken by extending the integration to the lower half plane 

as 

Tk = — f sign{y) z-'^dh , k>2, 

tTTK JC\D 

(3.2) 

Ti = - — / sign(?/) z'^d^z , 

t7C Jd 

where y = Im z. Contour integral representations are easily obtained using the Stokes 
theorem. They read 

Tk = \xm I yz-^ dz= \ i \y\z'^dz, k>l. (3.3) 
nk J-y mk JdD 

The non-closed integration contour 7 (shown in Fig. Hj) is the part of the boundary of B 
lying in the upper half plane (with the orientation from right to left). 

It is convenient to introduce the generating function of the moments T^: 

1 r d7' °° 
MM = -f ^^ = E^^^^'"S 1^1-^0 (3.4) 



■ni JdD z' ^ 



(here y' = Xmz'). The integral of Cauchy type in (13. 4p defines an analytic function 
everywhere inside D. In a small enough neighborhood of the origin this function is 
represented by the (convergent) Taylor series standing in the r.h.s. of (13.41) . 

For example, let B be the half-disk of radius R: \z\ < R, Xmz > 0. Then an easy 
calculation gives 



and 



M+iz) = -[2 + —- log- 



TT V Rz R + z 
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3.2 An electrostatic interpretation 



Similar to the standard harmonic moments from papers [2], [13], the moments have 
a clear 2D electrostatic interpretation. Let the interior of the domain B be filled by an 
electric charge with uniform density —1 and let the lower half plane (or just the real axis) 
be a grounded conductor. By the reflection principle, the electric potential ^~{z) in the 
upper half plane is equal to the potential created by the charge in B and the fictitious 
"mirror" charge of opposite sign in B: 



/ log 


z-z' 


Jb 


z — z' 



d'z'. (3.5) 



Let us show that the T^'s are coefficients in the multipole expansion of $ {z) in the 
interior of B near the origin. We have, for 2; G B: 

^ , d'^z' 1 r d'^z' 



1 f z'dz' 1 r z'dz' 

-z + 



2TTi JdB z' — Z 2Tli JdB z' — Z 

If y'dz' 1 f y'dz' 
z — z 



where we substituted z' = z' — 2iy' in the second hne and integrated the analytic parts 
by taking residues. Since the function under the integrals vanishes on the parts of the 
contours along the real axis, we can eliminate them and combine the two integrals into 
a single integral over 9(B U B) = dD: 

1 r \v'\dz' 

d,^~iz)=z-z - -i ^-^^^ = z-z-tY kTkz''-\ (3.6) 

71 JdD Z' — Z ^ 

where the second equality follows from (13. 4p . Since $~(x) = at real x, we obtain the 
expansion of the ^~{z) around in the upper half plane: 

^-{z) = hz--zr-tj:uz'-z'). (3.7) 

^ k>l 

Similarly, expanding ^'{z) around 00 in the upper half plane, we get: 

^,<l>-{z) = -^YVkZ-''-\ ^-{z)=iY^{z~^--z-'^), (3.8) 

k>l k>l 

where 

Vk = -Im I z^d^z (3.9) 
vr JB 

are moments of the interior. Their generating function, M_(z), is given by the same 
Cauchy-type integral (13. 4p for z outside D: 

1 r \n'\ dz' °° 
M_{z) = -l \yA^ = Y^VkZ-'-\ \z\^oo. (3.10) 
m JdD z' — z ^ 
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Sometimes an equivalent electrostatic interpretation appears to be more convenient. 
Let us assume that there is no conductor but the domain B in the lower half plane is 
indeed filled with the "mirror" charge. In this case formulas f l3.6p . fl3.7l) and fl3.8p admit 
continuation to the lower half plane which is achieved by complex conjugation of both 
sides. For example, complex conjugation of (13. 6p yields d^^'iz) = z — z + i J2k>i kT^z'''^, 
z E B, which can be rewritten as dz^~{z) = z — z + i J2k>i kT^z^"^, z E B. 



3.3 Local uniqueness of a fat slit with given moments 

Here, we prove the local uniqueness of a fat slit with given moments. The deformations 
changing only one moment will be constructed in the next subsection. 

For the purpose of this section it is convenient to work with symmetric domains 
D = BU B rather than with fat slits themselves. Let D{t) be a one-parameter deformation 
in the class of symmetric domains such that D(0) = D and dtT^ = for all k = 1,2, . . .. 
We shall show that any such deformation is trivial: D{t) = D(0) (at least in a small 
neighborhood of t = 0). To see this, consider the t-derivative of the function M^{z). A 
simple calculation (see Appendix A) shows that 

9,M,W = 1/ ^ign(y>.,.(.-) 

m JdD z — z' 

where Vn{z) = 6n{z)/6t is the "velocity" of the normal displacement of the boundary 
at the point z. (If x{a,t), y{(J,t) is any parametrization of the contour, then f„ = 

^{daydtx — daxdty), where dl = \dz\ = {dxY + {dyY is the line element along the 
contour and f„ is positive when the contour moves to the right of the increasing a 
direction). The Cauchy-type integral in the r.h.s. defines analytic functions both inside 
and outside the contour. For z inside the contour, choosing a neighborhood of such 
that 1 2; I < \z'\ for all z' G 9D, we can expand M^{z) as in (13.40 and find that 

00 

dtM+{z) = Y.kdtTkz''-^ =Q 

k=l 

for all z in this neighborhood. By uniqueness of analytic continuation dtM^{z) = 
everywhere in D. According to the property of integrals of Cauchy type this means 
that the function dtM+{z) is analytic in C \ D and is given there by the Cauchy integral 
(13. lip . The boundary value of this function is sign (y) Vn{z)\dz\/dz with realvn{z) almost 
everywhere on the contour (in our case actually everywhere except maybe the two corner 
points on the real axis). Furthermore, deformations of symmetric domains preserving 
symmetry with respect to the real axis obey the condition 

f sign{y)vn{z)\dz\ = 

JdD 

(because Vn{z) = t>„(z)), which means, according to the Cauchy integral representation, 
that the function dtM^{z) has a zero at infinity of at least second order. Invoking the 
technique from the theory of boundary values of analytic functions [T7], one can 
prove (see Appendix B) that any analytic function with such properties in C \ D must 
be identically zero. Therefore, Vn = 0, i.e., the deformation is trivial. 
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3.4 Special deformations of fat slits 



In order to define deformations tfiat cfiange one of tlie moments keeping all other fixed 
we need the odd Green's function for symmetric domains introduced in Section 2.2. 

Fix a point a G H \ B and consider small deformations 5~ of a fat slit B defined by 
the infinitesimal normal displacement of the boundary as follows: 

«0 = -|5nG-(a,0, ee7, aGe\B. (3.12) 

(here ^ G and e 0). These are analogs of the deformations 5an{^) = — |5nG(a, ^) 
from [13] which generate dToda flows in the space of compact domains in the plane. 
As we shall soon see, 6~ generate, in the same sense, dKP-fiows in the space of fat slits. 
Extending the definition of the 6~ to the lower half plane, we can define the corresponding 
deformation of the symmetric domain D = B U B: 

6-n{0 = -^sign{ImOdnG~{a,0, ^ e dD . (3.13) 

Clearly, 6~n{C,)/e as e ^ is to be understood as normal velocity of the boundary under 
the deformation. 

An important comment is in order. For the deformations 6~ to be well defined around 
the points x_, x+, we assume that the angles a^. between the curve 7 and the real axis 
are strictly acute: < a-t < 7r/2 (see Fig. [2]). Since p'{z) ~ (z — x±)"± around the 
corner points (for a rigorous proof see, e.g., Lemma 2.8]), it is seen from 02.81) that 
the normal velocity of the boundary near the corner points tends to zero a.s C, x± and, 
moreover, so does the angular velocity of the parts of the boundary near the corners 
(which is of order |p"(z)|). This means that the points x± and the angles a± remain 
fixed. For not strictly acute angles a+, «_ the situation is much more complicated. For 
example, the angles can immediately jump to other values and the deformations are not 
always well defined (cf. [H]). This case deserves further investigation. 

Expanding the Green's function as in fl2.10p . one can introduce the deformations 

S^'^'^n (0 = ^ sign (JmO dnImBkipiO) , ^ e dD . (3.14) 

Like S~, the deformations S^'^"^ do not shift the endpoints of 7. 

It is not difficult to show that changes the harmonic moment keeping all 

other fixed. Indeed, assuming that a G D, we write: 

5(^'-)M^(a) = 1 / sign(y)^(^.-)n(.) . d^JrnBM^ 

ni JdD a — z 2Tci Jdo a — z 

where the extra sign in the definition of S^^'~^ cancels the sign(?/) in the definition of 
M_|_(a). Because ImBk{p{z)) = on 5D, we have: 

27ri JdD z — a 2ni Jao z — a 

where we have used the fact that the function B^. — z^ is analytic in C \ D and vanishes 
at 00, and so does not contribute to the integral. We thus see that 6^^~''Tj = eSkj. 
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3.5 Vector fields in the space of fat slits 



Deformations which depend on 7 in a smooth way can be represented by vector fields in 
the space of fat shts. Let 6n{z) be any small deformation of a fat slit. Given a functional 
X on the space of fat slits, its variation reads: 

6X 



The variational derivative 6X/6n{C,) has the following meaning: 

SX , (5(^)X 
lim 



Here 5^^^X is the variation of the functional under attaching a small bump of area e — > 
at the point ^ G 7 (a symmetric bump is assumed to be attached at the point ^). Let 
5n{^) = eg{$,), e — 0, then we define the vector field (the Lie derivative) V'-^-': 

6X 



V-(a)X^~j -ilLc)„G-(o.{)|d{|. (3.15) 



Applying these general formulas to the deformations 5„ (see fl3.12p ). we can write 
5~X = eV^(a)X, where the vector field V~(a) acts on functionals as follows: 

1 r 6X 

MO 

This equation gives an invariant definition of the vector field V~(a) independent of any 
choice of coordinates. According to the Dirichlet formula (12. 5p . the action of V~(a) 
provides the harmonic extension of the function TT6X/6n{^) from 7 to EI bounded at 
infinity and equal to on the rays of the real axis outside B. 

In the local coordinates T^, V~(a) is represented as an infinite linear combination 
of the vector fields d/dTk which can be thought of as partial derivatives. To find it 
explicitly, we calculate 

5-Tk = -\ f Im{z''')6-n{z)\dz\ = f Im{z-'')dnG-{a, z)\dz\ = 

irk J7 irk Jj k 

where the last equality follows from the Dirichlet formula (12. 5p for symmetric domains. 
Now, given a functional X on the space of fat slits, and assuming that X is a function 
of the moments Tk only, we write 

k>iOlk k>l^ 

SO V~(a) is given by 

V-(a) = -2 5:ixm(a-^')A. (3.16) 

fc>i k 

The vector fields V~(a) are "half-plane" analogs of the vector fields V(a) introduced 
in [14] via their action on functionals X in the space of all domains: 

V{a)X = -\[ ^drrG{a,Om (3.17) 

2 JdD dn{t,) 

(cf. (I3.15P ). Below we will show that the dKP hierarchy is related to the vector fields 
V^(a) in the same way as the dToda hierarchy is related to the V(a). 
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4 The dispersionless KP hierarchy 



4.1 Lax equations 

The dispersionless KP (dKP) hierarchy is encoded in the Hadamard formula (12. 6p . To 
see this, fix three points a, 6, c G H \ B and find 5~G~{a, h): 

5-G-{a,h) = - ^ / dnG-{a,z)dnG-{h,z)dnG-{c,z)\dz\ . (4.1) 
47r J7 

This formula shows that the quantity V^(a)G^(6, c) is symmetric with respect to all 
three arguments: 

V-{a)G-{h, c) = V-{b)G-{a, c) = V"(c)G-(a, b) . (4.2) 
Using the expansions (12.10p and f l3.16p . we get 

^Im{Bi{piz)) = ±-Im{BMz)) , k,l>l. 
oTk oTi 



Since Bk{p(z)) = Bf^(p{z)), we can easily separate holomorphic and antiholomorphic 
parts of this equality and rewrite it as a relation between functions of z only: 

dBi{p{z)) dBMz)) 



on dTi 

In particular, 



(4.3) 



dBkipjz)) ^ dp{z) 

m m ' ^ ' ' 

Treating p rather than z as an independent variable and passing to the inverse map, z{p), 
one can bring this equality to the form 

dzjp) _ dBk{p)dz{p) dBk{p)dz{p) _ rp . X XI 
on ~ dp dT, dT, dp -WP),z[p)i. ^4.5j 

Recalling that Bk{p) = (z'^(p))>o (see f l2.12p ). we recognize the standard Lax equations 
of the dKP hierarchy. 

So far we assumed that p does not belong to the segment (see Fig. [3]). This 

segment is a branch cut of the function z{p). On this cut we can write 

z{p ± iO) = x{p) ± i\y{j)) \ , p- < p < p+ . (4.6) 

Here \y{p) \ = y{p + iO) = —y{p — iO) > 0. Outside the cut, the real- valued function x{p) 
has the same expansion (12.20 as the z{p). On the cut, it is given by the principal value 
integral 

.(p) = p+ip.v.rMi;W. (4.7) 

TT Jp_ p — p 

Because Bk{p) is a polynomial with real coefficients, one sees from (14.50 that x{p) and 
\y{p)\ obey the same Lax equations: 

^ = {B,ip),xip)}, ^^ = {5.(p), |y(p)|}, pe[p-,p^]. (4.8) 
In the next subsection we show that 2|y(p)| is the Orlov-Shulman function. 
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4.2 The Orlov-Shulman function 



Consider the functions M^{z), M_{z) defined by the integrals of Cauchy type (13 .4^ . 
fl3.10p for z inside and outside the domain D = B U B respectively. They can be also 
represented by the Taylor series 

M+(z) = mz'-' , M_{z) = ^ VkZ-''-' 

k>l k>l 

which converge in some neighborhoods of and cxo respectively. The function M+(z) is 
analytic everywhere in D while M_{z) is analytic everywhere in C \ D (with zero of second 
order at oo). Moreover, for analytic arcs 7 both M^{z) and M_{z) can be analytically 
continued across the arcs 7 and 7 everywhere except their endpoints on the real axis, 
where both functions have a singularity. Therefore, the function 

M{z) := M+{z) - M.{z) = ^ kT^z^'^ - ^ V^^-^-^ (4.9) 

fc>l k>l 

is analytic in a neighborhood of the boundary of D (excluding the points x± on the real 
axis) and, by the property of the Cauchy-type integrals, is equal to 2\Imz\ = 2\y\ on 
7 U 7. This can be also seen from formulas (13.61) . (13.81) (together with their extensions to 
the lower half plane) taking into account that the derivatives of the electrostatic potential 
are continuous at the boundary: 



out 



where the upper (lower) sign is taken for z in the upper (lower) half plane. We see that 
for z in the upper half plane M{z) is the analytic continuation of the function 2y from 
the contour 7, while for z in the lower half plane M{z) is the analytic continuation of 
the function —2y from the complex conjugate contour 7, i.e., M{z{p)) = 2\y{p)\. 

Let S{z) be the Schwarz function of the contour 7, i.e., an analytic function such that 
S{z) = ^ for 2; on 7 (see [20] for details). For analytic contours, it is known to be well 
defined in some strip-like neighborhood of the curve. Clearly, the Schwarz of the complex 
conjugate contour 7 is then S{z) = S{z). By uniqueness of analytic continuation, we can 
express M{z) in terms of the Schwarz function: 

i{S{z) — z), Imz > 
M{z) = { _ (4.10) 
—i{S{z) — z), lmz<0. 



Let us show that 



dT,M{z) = dA{p{z)). (4.11) 



Consider the change of the S{z) under the deformation . If 2; G 7, then, using the 

dxSiz) 

identity Vn(z) = , for the normal velocity of the boundary under a deformation 

2t^f¥h) 



with a parameter T, we can write 

6-S{z) = -2z,/s^)^d^G-{a,z). 
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Since \JS'{z) = \dz\/dz = 1/t{z), where t{z) is the unit tangent vector to the curve 7 
(represented as a complex number) and dnG~'{a,z) = —2iT{z)dzG~{a,z), we have 

6-S{z) = -2ed,G-{a,z) 

and so, 

\/-{a)S{z) = -2d,G^{a, z) (4.12) 

for z G 7 and, by analytic continuation, everywhere in the neighborhood where >S'(z) is a 
well defined analytic function. Expanding both sides as in (12.101) . (I3.16p . we finally get 

dT,S{z) = -%d,B^W)), (4.13) 

which is equivalent to (14. lip . 

An important particular case of (14.111) is 

dT,M{z) = d,p{z) . (4.14) 

Passing to partial derivatives at constant p, one can rewrite it in the form of the "string 
equation" : 

{z{p),2\y{p)\} = l, pe[p_,p+]. (4.15) 

This relation together with the Lax equations ( 14. Sp show that M{z{p)) = 2\y{p)\ is the 
Orlov-Shulman function [21] of the dKP hierarchy which describes deformations of fat 
slits. 

A closely related useful object is the indefinite integral of M{z). Using the notation 
of the previous subsection, we introduce the function 



n{z) = / M+{z)dz+ / M_{z)dz 

Jo Jz 

fe>l k>l 

It is analytic in the same strip-like neighborhood of the curve 7 where the Schwarz 
function is well defined. Since the electrostatic potential ^~{z) is continuous on 7, it 
follows from ([321), (M that 

Imn{z)=y^, ze-f. (4.17) 

The real part of Q{z) at z G 7 has the meaning of the partial area beneath the curve 7. 
More precisely, let 

nx 

A{z) = / y'dx' , z = X + iy G 

J X- 

be the partial area of B cut from the right by a line orthogonal to the real axis and 
passing through z, then 

nen{z{)~nen{z2)=2{A{z^)-A{z2)), ^1,22^ 7. (4.18) 

By construction, partial T^-derivatives of the function f2 at constant z are the generators 
of the fiows: 

Bu{p{z)) = dT,^{z) . (4.19) 

In this sense the function Vt = fl{z; {Tj}) solves the whole set of equations ( 14. 3p and thus 
provides a solution to the dKP hierarchy. 
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4.3 The string equation 

The string equation can be also derived in a different way along the lines of [I3]. The 
idea is to use equation (14. 2p with one of the three points lying on the contour 7: 

V-{OG-{a,b) = V-{a)G'{b,0, ^€7, 

and the other two points tending to infinity. From fl3.15p we see that the l.h.s. is equal 
to ^dnG^{a,^)dnG^{b,^) which is 2Im{a~^)Xm{h~^){dn1rnp{S,)Y as a,b 00. The 
r.h.s. in the same limit is —A1m{a~^)1m{h~^)dTiXmp{^). Equating them, we obtain the 
important relation 

2dT,Imp{z) = -\d,p{z)\'^, z e . (4.20) 
Its extension to the lower half of the boundary of D reads 

2dTjmp{z) = |9^p(z)p, z G 7. (4.21) 

Passing to the variable p with the help of the identity 

dT,p{z) 
dzP{z) 

we rewrite equations (14.201) . (I4.2ip in the form 

dz{p — ie) dz{p + ie) dz{p — ie) dz{p + ie) 



dp dTi dTi dp 



zsigne, e^O. (4.22) 



Note that the equation at e < is obtained from the one at e > by complex conjugation. 
Passing to the functions x{p) and y{p), we see that (14.220 coincides with (14.150 . 

Eq. (14.220 can be cast into the form of an evolution equation for z{p). To derive it, 
we start from the Hadamard formula for the deformation 6^^~^: 

6^'^-^G-{a,z) = ^ f dnG-{a,OdnG-{z,OdnImp{Om. 

Expanding the Green's function at \a\ 00 (see (12.101) ) and using the fact that 6''^'^^Ti = 
e, 5(i'-)Tfe = at > 2, we rewrite it as 



drjmpiz) = j-J dr.G~{z,0 WiOl'm 



which, after extracting the holomorphic part and passing to the integration in the ]9-plane, 
reads 

1 dp 



{p{z) -p)\z'{p)\ 



2 



In terms of the function z{p) we have: 



which is the desired evolution equation equivalent to (I4.22p . The latter is obtained from 
(I4.23P by taking the jump of both sides across the segment Because the function 

z{p) is analytic in the upper half plane, this is equivalent to the full equation (14.230 . 
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5 The dKP hierarchy in the Hirota form 



Here we reformulate the dKP hierarchy in the Hirota form using the dispersionless "tau- 
function" (free energy) and clarify the geometric meaning of the latter. 



5.1 Functional F and its variations 

Given a domain D, not necessarily symmetric, one can introduce the functional F: 

F = --J I \og\z-^-C'\d'zd\. (5.1) 
vr^ JD JD 

If the boundary is analytic, F is the "tau-function for analytic curves" introduced in [12] 
(more precisely, a properly rescaled logarithm of the dToda tau-function). In the 2D 
electrostatic interpretation, it has the meaning of the electrostatic energy of a uniformly 
charged domain D with a compensating point-like charge at the origin. It was found in 
that the Green's function G{a, b) is given by 



G{a, b) = log \a-^ - r^l + ^V(a)V(6)F, 



(5.2) 



where the vector field V(-2) in the space of all domains is defined in (13.17P (see [131 E 
for more details). We are going to derive an analog of equation (15. 2p for fat slits in the 
upper half plane. 

Given a fat slit B, consider the following functional: 



TT' JB JB 



log 



z-C 



(5.3) 



It has the meaning of the 2D electrostatic energy of the uniformly charged fat slit in the 
presence of a conductor placed along the real axis. Taking a variation of F~, it is easy to 
find how the vector field V~(a) acts on F~ . We use the general relations given in section 
4.2. We have: 



TT 



5F- 



- / log 

TT JB 



The function in the r.h.s. is already harmonic and bounded in EI \ B as it stands, hence 

2 





a — z 


/ log 


IB 


a — z 



d^z = $ (a) 



(5.4) 



Expanding both sides as a — * oo and comparing coefficients in front of the basis harmonic 
functions Im{a~^), we find: 



dF^ 



-Im / z'^d'^z 

TC JB 



(5.5) 



where Vk are the interior harmonic moments. Proceeding in a similar way, we find 



TC ■ 



6<l>-{a) 



2 log 



C, — a 
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The r.h.s. is harmonic (in ^) everywhere in EI \ B except the point a. This singularity 
can be canceled by adding the Green's function G~ (which vanishes on the boundary). 
Therefore, 

V-(fe)$-(a) = -21og ''"^ 



a — b 



+ 2G'{a, b) 



and we obtain the formula for G {a, b), 



G (a, b) = log 
which is a "half-plane" analog of fl5.2p 



a — b 



a — b 



-V"(a)V-(6)F- 



(5.6) 



5.2 Hirota equations for the dKP hierarchy 



Equation (15. 2p is known to encode the dToda hierarchy in the Hirota form. Equation 
(15. 6p does the same for the dKP hierarchy. To see this, let us apply the arguments from 

Combining (15. 6p and (12.71) . we obtain the relation 



log 



p{z) 


-p{z') 


2 

= log 


z — z' 


p{z) 




z — z' 



+ V-{z)V-{z')F- 



(5.7) 



which implies an infinite hierarchy of differential equations for the function F . Recall 
that the conformal map p{z) is normalized as 



u 



p(z) = z^ h0(l/2^), z ^ oo 



(5. 



(see (12.11) ). Tending 2;' — > oo in (15.71) . one gets: 

Xmp{z) = Xm -2+2 ^Ti'^ {z)F~ . 
The limit 2; — > cxd of this equality yields a simple formula for the capacity: 

u = dl^F- . 



(5.9) 



(5.10) 



Let us separate holomorphic parts of these equations, introducing the holomorphic part 
of the operator \/~{z): 



D{z) = Y.—dn, V-{z)=i[D{z)-D{z)]. 



k>l 



(5.11) 



Equation (15.70 then implies the relation 

p{z)-p{z') z-z' 
log -7-^ -7^ = log 



p{z) ~p{z') 



iDiz)V-{z')F- 



(5.12) 



which is holomorphic in z. In the limit 2;' ^ 00 it gives the formula for the conformal 

map p{z)\ 

p{z) = z + dT,D{z)F- (5.13) 
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(this formula also follows from fl5.9p ). In a similar way, equation (15.121) implies the 
relation 

z — z 

which is holomorphic in both z and z' . Taking into account (I5.13p . we rewrite it as 
follows: 

1 _ ^-Diam,)F- ^ _ m^md^^F^. (5.15) 

a — 

It is the dKP hierarchy in the Hirota form. We see that the function F~ is the disper- 
sionless tau-function for this hierarchy. The double integral representation (15.31) clarifies 
its geometric meaning. 



6 A growth model associated with dKP hierarchy 

The special deformations from section 4.1 suggest to introduce a growth model which 
is associated with the dKP hierarchy in the same way as the Laplacian growth [HI [9] 
of compact planar domains at zero surface tension is associated [5] with the dToda 
hierarchy. In fact the model to be introduced is also of the Laplacian type, i.e., the 
interface dynamics is governed by the Darcy law, but differs from the standard one by 
boundary conditions. 

The idea should be already clear from section 4.1: to consider growth of a fat slit 
under the deformation which changes only the first harmonic moment Ti keeping 

all other fixed and to identify Ti with time T. 

The corresponding growth problem can be formulated as follows. Consider a fat slit 
B(T) with a moving boundary 7(T), where T is time, and suppose that the motion of 
the boundary follows the Darcy law: 

Vn{0 = \dnm. i^l- (6.1) 

Here f„(0 = ^^iO/^'^ is the normal velocity of the boundary at the point C, and (p{z) is 
a harmonic function in H \ B such that 

(i) = on 7 and on the rays of the real axis [— oo, [xj^, +oo]; 

(ii) (j){z) = Im z + o(l) as Im z —>■ +oo. 

Clearly, (j){z) = Imp{z), where p{z) is the conformal map (12.11) . is harmonic in H[\ B and 
obeys these conditions. Comparing with (I3.14p at = 1, we see that the dynamics is given 
by the deformation 5^^'~^ at any point in time T = Ti and all the higher moments are 
integrals of motion. In other words, for our growth process 9tM+(z) = 1. Equivalently, 
the dynamics can be reformulated in terms of the inverse conformal map z{j), T) as the 
"string equation" in the form (14.221) . 



Xm 



dXp-tO,T)dTzip + iO,T)]=^, pe[p,,p+], (6.2) 
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Figure 6: The curve + = in the wpper half plane. The tangent lines at the 
origin are at an angle of A'b^ to the real axis. 

or in the form of the evolution equation fl4.23p (a similar equation for the Laplacian 
growth in the standard setting is well known [22]). As was already mentioned, the 
growth process is well defined if both angles a± between 7 and the real axis are acute. 
Then these angles and the points x_, x+ stay fixed all the time. 

Comparing this setting with the standard Laplacian growth in the upper half plane, 
we see that the conditions on are very similar if not the same: = on an infinite 
contour from left to right infinity, harmonic above it and tends to 1m z as Tmz +00. 
However, in our case, unlike in the standard one, only a finite part of the boundary 
(namely, the part which lies above the real axis) moves according to the Darcy law 
while the remaining part (the rays of the real axis) is kept fixed despite the fact that 
the gradient of is nonzero there (Fig. [5]). We do not know a proper hydrodynamic 
realization of this growth process. 

So far we assumed that x_ was strictly less than x+, so that the base of a fat slit 
was a segment of nonzero length. The setting of this section allows us to consider the 
degenerate case x_ = x+ = when the base of a fat slit consists of one point. The first 
harmonic moment Ti = T as well as the function Mj^{z) are still well-defined but Mj^{z) 
is singular at z = and thus can not be expanded into the Taylor series around this 
point (this means that the higher harmonic moments (13. ip are ill-defined). 

The simplest explicit solution to equation (16.20 known to us describes self-similar 
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growth of a "fat slit" with degenerate base. The function 

performs a conformal map from the exterior in HI of the curve |2;^ + T^| = T^, or, in polar 
coordinates, 

R{e) = rV-2cos2e, ^<e<^, (6.4) 

to the upper half plane. This curve is shown in Fig. O In this case x_ = = 0, 
a„ = a+ = 7r/4, p± = ±2T. The inverse map has the form 

z{p,T) = ^{j? - ^T^fl\v^{j? - AT^y/^y^'. (6.5) 

One can check that it does solve equation (16. 2p . More results on explicit solutions to 
Laplacian growth of fat slits will be published elsewhere. 



7 Concluding remarks 

We have constructed a parametric family of conformal maps of the upper half plane which 
is related to the dKP hierarchy with real "times" in the same way as conformal maps 
of the unit disk onto compact domains in the plane with smooth boundary are related to 
the dToda hierarchy with complex conjugate "times" tk, i^. Like in the dToda case, the 
deformations of domains ( "fat slits" ) in the upper half plane induced by dKP flows have 
a physical interpretation as Laplacian growth with certain type of sources or sinks at 
infinity. At the same time, our construction extends the well known connection between 
the dKP hierarchy and conformal maps of slit domains. In all cases, the conformal map 
plays the role of the Lax function. 

However, a limiting procedure from fat slits to usual slits is singular and is not easy 
to trace on the level of the Lax equations. We hope that a better understanding of this 
limit will further clarify the geometric meaning of solutions to equations of the dKP 
hierarchy. We also expect that yet more general solutions can be obtained by the same 
method applied to the case of a background charge distributed in the upper half plane 
with a non-uniform density, in accordance with a similar construction given in [23] . 

The solutions discussed in this paper have a nice geometric meaning but it seems to 
be very hard to express them analytically in a closed form. In this respect the situation is 
less favourable than in the dToda case, where some simple explicit solutions for conformal 
maps as functions of a finite number of nonzero harmonic moments (corresponding, for 
example, to a parametric family of ellipses) are available. It is clear that in our case the 
situation when only a finite number of the moments are nonzero can not be realized 
because the local behavior of their generating function M^{z) near the points x_, a:_|_ 
(which can be found from the integral representation fl3.4p ) shows that they are branch 
points of this function. This suggests that the corresponding solutions may be similar to 
the multi-cut solutions to Laplacian growth discussed recently in [21]. 
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8 Appendices 



Appendix A 

Here we give some details of the derivation of the formula (13.111) for the time derivative 
of the function M+{z): 

mM = ^l ^^g^(^)^"(^) \d,\ . (8.1) 

Til JdD a — z 

Here Vn{z) is the velocity of the normal displacement of the boundary at the point z. Let 
z{a, t) = x{a, t) + iy{(J, t) be any parametrization of the contour such that a is a steadily 
increasing function of the arc length, then it is a simple kinematical fact that 

Vn = ^{davdtx - daxdty) , (8.2) 



where dl = \dz\ = y (dx)^ + {dyy is the line element along the contour. According to our 
convention, Vn{z) is positive when the contour, in a neighborhood of the point z, moves 
to the right of the increasing a direction. 

Without loss of generality, we assume that cr varies from to 2tt, and, furthermore, 
for symmetric contours we choose a in such a way that z{2tt — a,t) = z{(J,t), z{0,t) = 
z{2Ti,t) = x+, z{n,t) = x_ and y{<J,t) is positive for < a < tt. We have: M+(a) = 
I (a) + /(a), where 

j(a) = }_ r y(^) ^^(^) 



wi Jo z[a) — a 

A straightforward calculation gives: 

■ rsj.. r ( dtyd„z + ydtd^z ydtzd„z\ 

iTfdtlia) = I da 

\ z — a z — ar / 



z — a 



r r-^y^^^±y^da+yd,zd(- ' 

Jo \ z — a \z 

r dtyd^z - djzd^y ^ y{a)dtz{a 
Jo z — a 



z[a] — a 



The last term obviously vanishes and we obtain, using (18. 2p : 

1 r Vn{z) \dz\ 



dJia) = - I 



7 a — z 



Adding dtl{a), we get (lO) . 



Appendix B 

In this Appendix we outline the proof of the following proposition. 
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Proposition 1. Let D be a compact domain bounded by a closed piecewise analytic 
contour F = dD in the plane with a finite number of corner points. Consider the 
function h{z) defined by the Cauchy-type integral 

ft(,) = _L/f«M!, (8.3) 



27ri Jr z — ^ 

where p(^) is a bounded real- valued piecewise continuous function on F such that 

£p(OMel = (8.4) 
and assume that h{z) = for all z E D. Then p = 0. 

One can try to prove this statement by means of the following elementary argument. 
Let r(^) = dC,/\d^\ be the unit tangential vector to the curve F at the point ^ represented 
as a complex number. If h{z) = for all 2; G D, then the properties of Cauchy-type 
integrals imply that piO/'^iO the boundary value of a holomorphic function h{z) in 
C\D vanishing at infinity. In fact h{z) is given by the same integral (18.31) . where z G C\D. 
Condition (18. 4p tells us that the zero at infinity is of at least second order. Let w{z) be 
the conformal map from C \ D onto the unit disk such that w(oo) = and r = lim zw(z) 
is real positive. By the well known property of conformal maps we have 



along the curve F. Therefore, 



and we thus see that 



\dz\ \dw\ 



r{z)=^\w'{z)\^^, (8.5) 
p{z) w'{z) 



is the boundary value of the holomorphic function h{z). Since w'{z) 7^ in C \ D, the 
function 

is holomorphic there with the purely imaginary boundary value 

P(^) 
i\w'{z)\ 

Then the real part of this function is harmonic and bounded in C \ D and is equal to 
on the boundary. By uniqueness of a solution to the Dirichlet boundary value problem, 
lZeg{z) must be equal to identically. Therefore, g{z) takes purely imaginary values 
everywhere in C \ D and so is a constant. By virtute of condition (18. 4p this constant must 
be which means that p = 0. 

However, this argument is directly applicable only for purely analytic contours for 
which all singularities and zeros of the function w'{z) lie strictly inside it. For contours 
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with corners, the corner points are singular points of the conformal map w{z). Some 
more work is required to make the above argument rigorous. Below we present another 
proof of Proposition 1, which makes use of some non-trivial facts about boundary values 
of analytic functions and actually works in a much more general setting than just a finite 
number of corner point^. It takes advantage of translating the proposition to a statement 
about analytic functions in the unit disk. 



Sketch of proof of Proposition 1. If h{z) defined by (18.31) is identically in D, then 
it is analytic in C \ D, is 0(1/ z'^) near infinity (because of (18.41) ) and has the boundary 
value p{z)/t{z) almost everywhere on F. (In our situation "almost everywhere" means 
everywhere except a finite number of points.) It is known [16, chapter III] that h{z) 
belongs to the Smirnov class E^{C \ D). 

Let z = ip{w) be the conformal map from the unit disk onto C \ D such that v^(0) = oo 
and ip{w) = r/w + 0(1) as w with real r. (The function (p{w) is inverse to w{z) 
introduced above.) Set h{w) = h{ip{w)). Since 

, , dz dipiw) ip'(w)dw '{''{w) 

r{ip{w)) = — — 7 = — — - — r- = — - — -— — - = iw 



\dz\ \d({){w)\ \(^'{w)\\dw\ \({)'{w)\ 
it follows from the above that h{w) is analytic in the unit disk with the boundary value 

Kw) = ^^^^^^^'^^i^^' (8.6) 
iw Lp'[w) 

almost everywhere on the unit circle and has zero of at least second order at w = 0. 
Clearly, the function H{w) = h{w)ip'{w) is analytic in the unit disk because the second 
order pole of ^p'iw) at w = is canceled by the zero of h{w). Furthermore, according to 
the Keldysh-Lavrentiev theorem [171 chapter 10], H{w) belongs to the Hardy class H^. 
But then the function H{w) = wH{w) belongs to the same Hardy class and takes purely 
imaginary boundary values 

-ip{<^{w))\<^'{w)\ 

almost everywhere on the unit circle. The characteristic property of functions from the 
class is that they can be represented by the Poisson integral of their boundary values 
with real positive Poisson kernel (see [161 chapter II, §5] or Theorem 3.9 in [17]). This 
means that H{w) must be purely imaginary everywhere inside the unit disk and hence 
must be a constant. Since H{0) = 0, the constant is 0, so h{z) vanishes identically and 

The proof extends word for word to a more general case when p is an integrable 
function with respect to \dz\ (not necessarily bounded) and the boundary of D is a 
rectifiable Jordan curve. It is crucial that the differential p{z)\dz\ is real valued on the 
boundary. If one dropped that assumption, the statement is false. Moreover, functions 
from the class E'^ can have real boundary values in domains with cusps (see an example 
in [25]). 



thank D.Khavinson who suggested this proof and explained it to me. 
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